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FIG. 12 is a graph representing preferred values of prefixes ment of the invention, the weighting factor a^. differs from one 

used in the system of FIG. 1. symbol to another, the elements of a given vector being 

multiplied by the same weighting factor. With an OFDM 

DETAILED DESCRIPTION OF THE PREFERRED modulator in the transmitter fimctioning in this way, blind 

EMBODIMENT 5 channel estimation in the receiver can be done simply and at 

low arithmetical complexity. In particular, the receiver can 

FIG. 1 shows an OFDM communication system in accor- constantly estimate and track the channel impulse response 

dance with one embodiment of the invention comprising a without any loss of data bandwidth. Moreover, the demodu- 

transmitter comprising an OFDM modulator 1 and a receiver lator at the receiver can have advantageous characteristics, 

comprising an OFDM demodulator 2, the transmitter and the lo ranging from very low arithmetical cost (at medium perfor- 

receiver coromunicating over a communication chamiel 3 . mance) to high arithmetical cost (very good system perfor- 

An input bit-stream b„€(0, l),n=0, 1 , . . . ,K-1 is modulated mance). 

onto a set ofN carriers whose carrier amplitudes are given by More particularly, in the preferred embodiment of the 

the vector X(k)=(Xo(k),Xi(k), . . . ,Xjv_i(k))^, corresponding invention, the prefix of D samples that is added in the guard 

to OFDM symbol number k. Afterwards, the time domain 15 interval comprises a pre-calculated suitable vector P^= 

OFDM signal is generated by means 4 which performs an (cq, . . . ,c^_i)^ of D samples that is independent of the data 

Inverse Fourier Transform operation, or preferably an Inverse and that is weighted by a pseudo-random factor a^^ that only 

Fast Fourier Transform ('IFFT') operation [F^]"^=[F^]^with depends on the number k of the latest OFDM symbol: 

[F^]^=([F^]0* where (•) ^ is the transposition operator and .^s^fj^.-f^^ „^ ^ .^w (k) 

(•)* is the complex conjugate operator: 20 {k)f ' ' ' " ' " '^^"^ Equation 2 

For the purposes of the analysis below, a second prefix/ 

x{k) = [F^r^Xik) = (xoik), xiik), ... , XN-i{k))'^., Equation 1 OFDM symbol vector is defined as follows: 

where x^^''"^'^''''\k)={xo{k)Mk), . . . ^M-i(k),Ci,,, 

^ 25 . . . ,Oik+iCn-L)^- Equation 3 

'■^^ •' ^ ^ ^o^i^N-i,o^^N-i ' Several choices for a^^ are possible. It is possible to choose 

^ a^^C, that is to say that a^^ can be of any complex value, 

and Wn = e'^'w . However, any a^^ with la^^M leads to performance degrada- 

3Q tion compared to preferred embodiments of the invention. 

_ , . 1, i . i X . -, It is possible to limit the choice of a^, somewhat less 

The resulting parallel signal x (k) vector is converted to a generally to a^eC with la^l=l . This choice usually leads to 

senes signal by a parallel-to-senes converter 5 a prefix, rep- ^ performance, but the decoding process risks to 

resented by the Dxl vector P^= c„, . . . ,c^_ J being inserted unnecessarily complex. 

mtothesisnalas guard mterval between each Of DM symbol ^ „ . , 

. 1 .° . / . 35 Accordingly, m the preferred embodiment of the present 

to produce a senes digital signal x„. ine senes digital signal . . ,° r- . i 

. .1 . 1 • 1 u\ u •+ 1 + mvention, the phase of is chosen so that 

x„ is then converted to an analogue signal x(t) by a digital-to- ^ t- k 

analogue converter 6 and transmitted over the channel 3. 

The channel 3 has a Channel Impulse Response H(k)=C(k) ^ 

and also introduces noise v. of* = e^'^ '", 

40 

At the receiver 2, an analogue signal r(t) is received and 

converted to a digital signal r„ by an anal ogue-to -digital con- 
verter 7. The digital signal r^ is then converted to a parallel where m is an integer, N is the useful OFDM symbol size and 
signal by a series-to-parallel converter r(k) and equalised and D is the size of the pseudo-random prefix. This choice is 
demodulated by equalisation and demodulation means 9 to particularly advantageous when using the specific decoding 
produce demodulated signals s''^^^(k). In the following analy- methods described below. 

sis, consideration of noise is omitted for the sake for simplic- For the sake of simplicity, the following analysis assumes 

ity. However, including the consideration of noise does not that the weighting factor has been chosen as 

significantly modify the results. 

In some known OFDM communication systems, the guard 

interval is used to add some redmidancy (D samples of redun- ^ _ ^.^.m 

dancy are added) by introducing a cyclic prefix, for example ak - e + 
in the following manner: 

x^^\k)=(x^_2){k), . . . :pcN-i(k):pco{k)yXi{k), . . . jji Integer. However, it will be appreciated that the mathemati- 

^^^'^ ■ cal adaptation to any of the cases presented above is straight- 

In other words, data from the end of the frame is repeated by ' j—i 

the transmitter in the guard interval to produce a prefix. proves to be very useful to choosdajsuch that its phase 

J , '.u ^u- I. A' ^ f^u +• changes from OFDM symbol to OFDMsymbol. The constant 

In accordance with this embodiment 01 the present inven- ^t^. n ^ ^ - • 

..1 .1 1 • «^ J A' ^ ^ f .r. prefix P^ IS preferably chosen with respect to-certamcnteria, 

tion, however, the prefix samples mserted as guard interval 01 60^ iir^n - 

OFDM symbol number k, a^.Co to a,.c^ „ are deterministic ^^^P^^ following: 

and are known to said receiver as well as to said transmitter. fi-equency domain, P^ is as flat as possible over the 

Theprefixes comprise a vector P^Kcq, • • • ,c^i)^of size Dxl frequency band used for data earners. 

that is common to the symbols multiplied by at least one In the frequency domain, P^ is as near to zero as possible 

weighting factor a^^, so that the prefixes have the overall form 65 for all unused parts of the band. 

a;fCo to a;j,.c^_i. The weighting factor a^^ may be constant In the time domain, P^ has a low peak-to -average-power- 

from one symbol to another. However, in a preferred embodi - ratio (PAPR) . 
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Thus, this method works very well if the prefix-spectrum is are shown in FIG. 12 by way of example, for the following 
non-zero everywhere in the FFT^^^ domain (and, of course, OFDM parameters: 

everywhere well above channel noise). This can be a trouble- Size of the Prefix in Time Domain: D=16 Samples 

some limitation in other circumstances. Size of the OFDM symbols in the frame: N=64 Samples 

A second embodiment of a method of channel impulse 5 Carriers where channel coefficients are to be estimated 
response estimation on D carriers in accordance with the (over N+D=80 carriers): Carriers 1 to 52 

present invention avoids this limitation, at the expense of Out-of-Band region: Carriers 76 to 80 
increased arithmetic cost. This second method does not esti- Maximum PAPR has not been limited 

mate fi^ based on a de-convolution in the FFT^,^ domain as Out-of-Band Radiation as low as possible 

presented above, but estimates FFT(^^^^^(^^^)((h^ 0^)) 10 Spectral Flatness as good as possible 
directly based on the received vector ([Ho]+[Hi])-Pr,. This is rr^i , i • j i i i .i 

•i 1 1 1 ... .1 1 . . The channel estimation is done by calculating the expec- 

possible by exploiting the observation: . , i ^ i pi . i 

tation value over a number oi samples oi the received vector 

[H(^N+D)x(N+D)\'iPD^^N^)^=(E/^o^^N-D^)^ Equation 13 as explained above. If the tracking of the chamiel is done 

. ^. . ^ J . J ^ -1 • T^T^ T ic based on a first estimation h(k-l)=rF;^l^- H(k-1) of the chan- 

This equationis represented in more detail m FIG. 10. In 15 ^ ^l x-rM-T^A/r i. i ^i. 

^1 . J .1 J ^1 1 1-1 • nel impulse response and a number R 01 OFDM symbols, the 

this second method, the channel impulse response is esti- . ^. . • j + j ^- n 

^ J . X- 11 . ^ ^ ^ first estimate IS then updated as follows: 

mated using the following steps: ^ 



Perfonn a FFT(^^n)v(Ar..7.) on 



Tf\ I\T 20 ^(^) = 



Perform a FFT^^^^^^^^^^^ on V^=(P^^O^^- 
Perform a component-by-component division frc- n x ^„ Jv^/rc-i x 

If desired, perform an IFFT on B-fj^.n^^' ^^r^rrr 

jppT m h{k) = [FofH{k). 

^ ^ (N+D)x(N+D)\"-N+D ) 25 

The last step of the list presented above is not essential for 

the basic equalization algorithm but may be useful, for ^ . ^ . . ^ . , . . 

example in algorithms used to reduce noise levels. ^^^^^ ""^^ ^^^^ °f the first method for channel estimation 

^ . .111 1 J -1 J -.1 X- . that has been presented above. Alternatively, the second 

I lie above methods have been descnbed with reierence to method can be a lied b 

the specific case where is constant and equal to 1. In 30 



preferred embodiments of the invention, however, the weight 

o.j^ of the prefix to each symbol k is a preferably complex 

pseudo-random factor that only depends on the number k of ^^^^ = 

the latest OFDM symbol. The adaptations to this method of ( r-\ 

the basic equations (shown in FIG. 9) are shown in FIG. 11.35 ■//(/:- 1) + \Fn^d\ ■ X ■ - 4 - n), ol_J 

It is found that equations 4 and 8 are to be adapted as 

follows: (S)([FN^D]-{PLONy 



hik) = [FN^D]"H{k). 



( {ak-PD)\ Equation 14 40 

To ^ I 



where the factors s„,n=0, 1 , . . . ,R-1 are positive real numbers 
" ^ ^ ^ ' ^ ■ that are used for normalization and weighting of the different 

/ {ak+i ■ Pd)\ Equation 15 contributions. Thus, for example it is possible to take older 

Ea,4 = E{[Hi]tM) + e[[Ho] ■ — ^— — J 45 OFDM symbols less into account for the channel estimation 



than later ones. The Fourier matrix [F] can be chosen in the 
- [Ho] ■ Pd- N+D carriers or D carriers domain [] 

Several equalization methods are advantageous using the 
The procedures for blind channel estimation described pseudo-random prefix OFDM. In general, the different meth- 
above remain applicable by setting Eo=E^o and E^=E^^,. 50 ods offer different performance-complexity trade-offs. 
This amounts to weighting the preceding and following' D A first embodiment of a method of equalization uses zero 

prefix-samples of each received symbol by the corresponding forcing in the N+D Domain and offers low complexity equal- 
^ or a^r.^ ^ ~ ^ respectively. ization. 
The values of the prefixes a^r^-P^ are chosen as a function of With 
selected criteria, as mentioned above. Values that have been 
found to give good results with the criteria: 

Low Peak-to-Average-Power-Ratio of the time domain A = 



60 



1 Of/t+l 

Signal 

Low Out-of-Band Radiation, that is to say maximise the 

energy of the prefix over the usefiil band and not waste Channd Impulse Response matrix can be represented as 

prefix energy over null carriers follows* 
Spectral Flatness, e.g. SNR of each channel estimates shall 
be approx. constant 

Low-Complexity Chamiel Estimation, i.e. by prefix spec- 65 ^ ^^^^^j + . ^^^^^^ ^ ^^^^^^ _^_ . ^^^^^^ Equation 16 

trum whose spectral contributions are mainly just phases 
(i.e. of constant modulus). 



11 
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-continued 



L ^DxN J L L ^DxN 



Thus, the equalized resulting vector is 



Equation 20 



[Ar={A"A)-U". 

where [A] is a rectangular matrix. 



Equation 21 



The invention claimed is: 

1. A method of communication using Orthogonal Fre- 
quency Division Multiplexing ('OFDM"), the method com- 
prising the steps of: 

generating bit streams b„ 6(0,1), n=0,l, . . . ,K-1 and the 
corresponding sets of frequency domain carrier ampli- 
tudes (Xo(k) to X^k)), where k is the OFDM symbol 
number, modulated as OFDM symbols to be transmitted 
from a transmitter, 

inserting prefixes as guard intervals in said sample streams, 

transmitting said OFDM symbols from said transmitter to 
a receiver, 

using information from said prefixes to estimate the Chan- 
nel Impulse Response (H^*^^) of the transmission chan- 
nels at the receiver, where (H^^^) is the length D vector 
defined as the channel impulse response vector in the 
frequency domain, denoted by superscripted, F, and 

using the estimated Channel Impulse Response (H^^^) to 
demodulate said bit streams in the signals received at 
said receiver, wherein said prefixes (a^ Cq to c^_i) 
are detenninistic and are known to said receiver as well 
as to said transmitter, where c is the set of vectors con- 
taining constant postfix samples, (a^) is a weighting 
factor proportional to 



55 



where[j|is the square root of -1 , N is the useful OFDM 
symbol size, D is the size of the prefix vector, and m is an 
integer, and fiirther performing the multiplication by a 60 
matrix proportional to 



R^^\k) = ^|N+D ■[V]-r{k), 
where 



[V] 



2n 
N+D 



10 



15 



where is the ratio of consecutive a^^ 



Pk = - — , 



The definition of the Moore-Penrose pseudo -inverse is, 
among others, discussed by Haykin in the book: "Adaptive 
Filter Theory" by Simon Haykin, 3'"'^ edition, Prentice Hall 
Information and System Science Series, 1996. Haykin uses 
the common definition 20 



calculating the frequency shifted CIR coefficients 



f.Skifted,F _ 



I N + D N + D ) 



25 



30 



35 



performing a component-by -component division 



i^'\k)=R^'\k)=is>H::T", 



performing a multiplication by a matrix proportional to 



I^^\k) = [V] 



^|N + D 



■I^^\k), 



40 



45 



50 



S(^^^(k)=[^-S^e 



extracting the N equalized samples corresponding to the 
k^^ data symbol to the vector S^^(k), and 
transforming the symbol s(k) int o freauencv domain bv 
performing a Fourier Transform 

(k). 

2. A method of communication as claimed in claim 1, 

wherein said prefixes (a^r^c^ to o-kCj^-i) comprise a vector ( 
P^) that is common to said symbols multiplied by at least one 
weighting factor (a;r,). 

3. A method of communication as claimed in claim 2, 
wherein said weighting factor (a^) differs from one symbol to 
another but the elements of a given vector (P ^) are multiplied 
by the same weighting factor. 

4. A method of communication as claimed in claim 3, 
wherein said weighting factor (a^) has a pseudo-random 
value. 

5. A method of communication as claimed in claim 1, 
wherein said weighting factor (a^^) is a complex value. 

6. A method of communication as claimed in claim 5, 
wherein the modulus of said weighting factor (a;^) is constant 
fi*om one symbol to another. 

7. A method of communication as claimed in claim 1, 
wherein estimating said Channel Impulse Response 



65 



comprises performing a Fourier Transform on a first 
vector (Vffj^) that comprises the received signal compo- 



